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Abstract. By recognizing the vital importance of two-hole Cooper pairs (CPs) in addition to the 
usual two-electron ones in a strongly-interacting many-electron system, the concept of CPs was 
re-examined with striking conclusions: namely, they are gapped and linearly-dispersive resonances 
with a finite lifetime — but provided the ideal-gas Fermi sea is replaced by a BCS-correlated unper- 
turbed ground-state "sea." Based on this, Bose-Einstein condensation (BEC) theory has been ge- 
neralized to include not boson-boson interactions (also neglected in BCS theory) but rather boson- 
fermion (BF) interaction vertices reminiscent of the Frohlich electron-phonon interaction in me- 
tals. Instead of phonons, the bosons in the generalized BEC (GBEC) theory are now both parti- 
cle and hole CPs. Unlike BCS theory, the GBEC model is not a mean-field theory restricted to 
weak-coupling as it can be diagonalized exactly. It reproduces the BCS condensation energy ex- 
actly for any coupling, and each kind of CP is responsible for only half the condensation energy. 
The GBEC theory reduces to all the old known statistical theories as special cases — including the 
so-called "BCS-Bose crossover" picture which in turn generalizes BCS theory by not assuming 
that the interelectronic chemical potential equals the Fermi energy. Indeed, a BCS condensate is 
precisely the weak-coupling limit of a GBE condensate with equal numbers of both types of CPs. 
With feasible Cooper/BCS model interelectonic interaction parameter values, and even without BF 
interactions, the GBEC theory yields transition temperatures [including room-temperature super- 
conductivity (RTSC)] substantially higher than the BCS ceiling of around 45K, without relying on 
non-phonon dynamics involving excitons, plasmons, magnons or otherwise purely-electronic mech- 
anisms. The results are expected to shed light on the experimental search for RTSC. 
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INTRODUCTION 

Boson-fermion (BF) models of superconductivity (SC) as a Bose-Einstein condensation 
(BEC) go back to the mid-1950's yj]-||4|], pre-dating even the BCS-Bogoliubov theory 
OD-lLZl]. Although BCS theory only contemplates the presence of "Cooper correlations" 
of single-particle states, BF models Illl-lfll.M- lfTill posit the existence of actual bosonic 



CPs. In spite of the central role played by CPs in both low- and high-r c superconduc- 
tivity, however, they are poorly understood. The fundamental drawback of early 
BF models, which took two-electron (2e) bosons as analogous to diatomic molecules in 
a classical atom-molecule gas mixture, is the notorious absence of an electron energy 
gap A(r). "Gapless" models are useful in locating transition temperatures if approached 
from above, i.e., T > T c . Even so, we are not aware of any calculations with the early 
BF models attempting to reproduce any empirical T c values. The gap first began to ap- 



pear in later BF models iHj-d. With two lflHllfl2Tl exceptions, however, all BF models 
neglect the effect of hole CPs included on an equal footing with electron CPs to give 
a GBEC theory consisting of both bosonic CP species coexisting with unpaired elec- 
trons in a ternary BF model. Although magnetic-flux-quantization measurements have 
established the presence of pair charge carriers in both conventional 1 151- lfl6ll as well as 
cuprate [fivh superconductors, no experiment has yet been done to our knowledge | fl~8^ | 
that distinguishes between electron and hole CPs. 

The "ordinary" CP problem [fl9|l for two distinct interfermion interactions (the 5-well 
Sllzil] or the Cooper/ BCS model lUllQH] interactions) neglects the effect of two-hole 
(2h) CPs treated on an equal footing with 2e-CPs — as Green's functions (220, on the 
other hand, can naturally ensure. However, a crucial confirmed result [|12ll is that the BCS 
condensate is a very particular BE condensate with equal numbers of 2e- and 2h-CPs, 
each contributing to one-half the condensation energy ||2~3ll . This was already evident, 
though not fully appreciated, from the perfect symmetry about £ = fx, the electron 
chemical potential, of the well-known Bogoliubov [|24jl v 2 (e) and u 2 (e) coefficients, 
where £ = h 2 k 2 /2m is the electron energy and m its effective mass. The GBEC theory 
(also appropriately viewed elsewhere as a "complete boson-fermion model," or CBFM) 
"unifies" [13] both BCS and ordinary BEC theories as special cases, and predicts 
substantially higher T c 's than BCS theory with the same Cooper/BCS model interaction 
that mimics electron-phonon dynamics. 



SIGNIFICANCE OF THE COOPER INSTABILITY 

A Bethe-Salpeter (BS) many-body equation (in the ladder approximation) treating both 
2e and 2h pairs on an equal footing reveals that, while the CP problem [based on an 
ideal Fermi gas (IFG) ground state (the usual "Fermi sea")] does not possess energy 
solutions with a real part, it does so when the IFG ground state is replaced by the BCS 
one. This is equivalent to starting from an unperturbed Hamiltonian that is the BCS 
ground state instead of the pure-kinetic-energy operator corresponding to the IFG. The 
remaining Hamiltonian terms are then assumed amenable to a perturbation treatment. 
As a result: i) CPs based not on the IFG-sea but on the BCS ground state survive 
through a nontrivial solution as "generalized" or "moving" CPs which are positive 
energy resonances with an imaginary energy term leading to finite-lifetime effects; 
ii) as in the "ordinary" CP problem, their dispersion relation in leading order in the 
total (or center-of-mass) momentum (CMM) hK = ft(ki +k2) is also linear (originally 
reported without proof in Ref. tf25|] . p.33)rather than the quadratic h 2 K 2 /2(2m) of a 
composite boson (e.g., a deuteron) of mass 2m moving not in the Fermi sea but in 
vacuum; and iii) this latter "moving CP" solution, though often confused with it, is 
physically distinct from another more common trivial solution sometimes called the 
Anderson-Bogoliubov-Higgs (ABH) |@](|0] p. 44), HQ collective excitation. The 



ABH mode is also linear in leading order and goes over into the IFG ordinary sound 



mode in zero coupling. All this occurs in ID 112811 . 2D H29I1 as well as in the 3D study 
outlined earlier in Ref. [|3Qh . In this section we focus on 2D because of its interest [31] 
for quasi-2D high-r c . cuprate superconductors. In general, the results will be crucial for 
BEC scenarios employing BF models of superconductivity, not only in exactly 2D as 



with the Berezinskii-Kosterlitz-Thouless (BKT) [32] [33] transition, but also down to 



(1 + e)D which characterize the quasi- ID organo-metallic (Bechgaard salt) [34]-[3 



and most recently multi-walled carbon nanotube [37], SCs. 

If K = ki + k.2 is the CMM and k = ^(kj — k2) the relative momentum wavevec- 
tors of the 2e bound state, and $k = E\ + £2 is its energy with E\ and £2 the 
energies of electrons 1 and 2, one uses the bare one-fermion Green's function 
Go(ki =K/2 + k,£'i = $k/2-\-E) for particle 1, and similarly for particle 2, where 
E = \{E\ —E%). The solution of the complete BS equation based on the IFG unper- 
turbed state with both 2e- and 2h-CPs included, and formed via the Cooper/BCS model 
interaction, is 

S Q = ±i2hco D / vVA - 1 ( 1 ) 

where X = VN(Ep) with N(Ep) the electronic density of states (DOS) for one spin, 
while V is a positive coupling constant and hato an energy cutoff, both defined below 
in (PT2)) . As the CP energy is pure-imaginary there is an obvious instability of the CP 
problem when both type pairs are included. This result was originally derived in Refs. 
no p. 44 and [38]; also, it was guessed in Ref. [12511 p. 167 without explicit mention of 



hole-pairing. It contrasts sharply with the familiar solution [11911 for 2e-CPs only, namely 



So = -2hco D /(e 2/X - 1) — ► -2h(0 D e- 2/X . (2) 

X — >0 

The first expression is exact in 2D and a very good approximation otherwise if hoop -C 
Ep, where (Oo is the Debye frequency. The sometimes misnamed "Cooper instability" 
© merely represents a negative-energy, stationary-state (i.e., infinite-lifetime) bound 
pair. We suggest, however, that unlike the apparent negative-but-real-^o "instability" © 
the genuine Cooper instability is really ([I]) so that the original CP picture is meaningless 
if2e- and 2h-CPs are treated on an equal footing, as consistency demands, since it leads 
to a purely-imaginary eigenvalue Sq. 

However, a BS treatment of pairs referred not to the IFG sea but to a BCS-correlated 
ground state "sea" vindicates the CP concept in terms of a new nontrivial solution. This 
is tantamount to starting not from the IFG unperturbed Hamiltonian but from the BCS 
one. Its physical justification is reinforced through the recovery of three expected results: 
a) the (trivial) ABH sound mode solution; b) the BCS T = gap equation; and c)finite- 
lifetime effects of a "moving-CP" nontrivial solution in either 2D [29] or 3D Q30D . Thus, 
the IFG Green function Go (ki , E\) is replaced by the BCS one, say, Go(ki ,E{) that now 
refers to an energy E\ = E^ = a/^ 1 2 + A 2 with ^ = h 2 k 2 /2m — Ep and A the T = 
fermionic gap. It also contains the Bogoliubov functions |[24ll vf. = |(1 — ^k/Ek) and 
u\=\ — v?. There are then two solutions to the BS equations. A trivial solution is the 
ABH energy eigenvalue $k, which when Taylor-expanded about K = gives for small 
X in 2D 

g K = ^K + 0(K 2 )+o(X), (3) 
V2 

where o{X) denote interfermion interaction correction terms that vanish as A — > 0. Note 
that the leading term is just the ordinary sound mode in an IFG with sound speed 
Vf I \fd in d dimensions, as determined straightforwardly from standard thermodynamic 



formulae. Secondly, there is a nontrivial moving CP solution of the BCS -correlated- sea- 
based BS treatment, which is entirely new and leads to the pair energy Sk which in 2D 
is [Hi 



±£ K = 2A+—hv F K+ 1 -^e 1 ^K 2 
2n 9 ko 



X 

—hv F K + 
K 12 k D 



1 hv FJ/X K 2 



+ 0{K 3 ). (4) 



Here, the upper and lower signs refer to 2e- and 2h-CPs, respectively, and 
ko = CQd/vf with COd the Debye frequency. A linear dispersion in leading or- 
der again appears, but now associated with the bosonic moving CP. From © 
the positive-energy 2p-CP resonance has an energy width Tk and a lifetime 

% K = h/2T K = h/2 \(X/n)hv F K+(hv F /12k D )e l / l K 2 that diverges at K = 0, falling 

to zero as K increases. Thus, "faster" moving CPs are shorter-lived and eventually break 
up, while "non-moving" K = ones are in infinite-lifetime stationary states. 



GENERALIZED BEC THEORY OF SUPERCONDUCTORS 

The GBEC theory [IT], [12] is described in d dimensions by the Hamiltonian H = 
Hq +Hi, lt . The unperturbed Hamiltonian Hq corresponds to a non-Fermi-liquid "normal" 
state which is an ideal (i.e., noninteracting) ternary gas mixture of unpaired fermions and 
both types of CPs, two-electron (2e) and two-hole (2h), namely 

H ° = s a khS + Y,E+(K)b+b K -Y,E-(K)c+c K (5) 

k u s K K 

where as before K = ki + 1*2 is the CP CMM wavevector while ££, = h 2 k 2 /2m, e.g., 
are the single-electron, and E±(K) the 2e-/2h-CP phenomenological, energies. Here 
a ks ( a k,j) are creation (annihilation) operators for fermions. Similarly (Z?k) and c£ 
(ck) are such for 2e- and 2h-CP bosons, respectively — although we do not attempt to 
construct them starting from the fermion operators. Two-hole CPs are postulated to be 
distinct and kinematically independent from 2e-CPs, all of which provides a ternary BF 
gas mixture. This postulate is firmly grounded on the experiments |l5ll- |[l7ll cited before. 
The operator Hq then has diagonal form and its exact eigenstates can be numerated by 
the sets of occupation numbers {...n^ s ...Nk---Mk...}. The occupation numbers n^ s of 
one-fermion states each take on only the two values and 1 , while those of the one-boson 
momentum-K states of 2p-CPs Afc, and of 2h-CPs Mk, take on all values 0, 1, 2, • ■ 
The exact eigenstates of the Hamiltonian Hq are then 



...n Ks ...N K ...M K ...) = nC^IT^r i b i) U^={ c k) I <>) (6) 




where | O) is the vacuum state for fermions and simultaneously for 2e-CP and 2h-CP 
creation and annihilation operators. Specifically, 

a Ks \O)=b K \O) = c K \O)=0. (7) 



If Nop is the number operator of the total number of electrons, charge conservation 
implies that 

Nop = £ a+ a kuSi +2 J>+£ K -2£c+c K . (8) 

k l5 S! K K 

If jl is their chemical potential, the exact eigenvalues of Hq — /lN op are 

E.., lks ...N K ...M K ... = [E+(0)-2n]N + [2n-E-(0)]Mo + Y,(e k -ti)n^ 

+ £ [E+(K)-2n]N K + £ [2n-E-(K)]M K . (9) 

K^O K/0 

The interaction Hamiltonian Him part of H = Hq + Hi nt consists of four distinct BF 
interaction vertices each with two-fermion/one-boson creation or annihilation operators, 
depicting how unpaired electrons (subindex +) [or holes (subindex — )] combine to form 
the 2e- (and 2h-) CPs assumed in the <i-dimensional system of size L, namely 

H * = L ^ /2 E/+W{< + >K,T a -k+^Ka Z,K + a -k + lK,I a k+lK,T^K} 

+ /- '''I/ C^K.Ik. " k.lK.. r K •« k lK./'k .K. <k} (10) 
k,K 

where k = i(kj — is the relative wavevector of a CP. The energy form factors 
f±(k) in (fTOl) are essentially the Fourier transforms of the 2e- and 2h-CP intrinsic 
wavefunctions, respectively, in the relative coordinate of the two fermions. The GBEC 
Hamiltonian H = Ho+Hi nt is very different from the well-known BCS Hamiltonian 

H BC =Hq +H? nt = £ E h al uS a h ^ + £ V kl ^a^a^a^a^ (11) 

ki,*i k^li 

with, say, the Cooper/BCS model interaction, with V > 0, 

f -V/L d if n-hG} D <e kl ,e h <n + h(o D n? . 
kl ' 1[ 1 otherwise. K ' 



In Refs. Illllllll2ll the energy form factors f±(k) for the GBEC theory are taken as 

f .(e) = Sf { f ^[E±(0) - 8e] < e < \[E ± (0) + 8e] 

J 1 ' \ otherwise v 

in order for it to reduce properly to BCS theory with (fT2)) in the limit to be explained 



below, provided one identifies / and 8e with ^J2Vhtt>u and hap, respectively. One 
then introduces the quantities Ef and 8e as new phenomeno logical dynamical energy 
parameters (in addition to the positive BF vertex coupling parameter /) that replace the 
previous such E±(Q), through the definitions 

E f = -[E+(0)+E-(0)] and 8e = -[£+(0) -£-(0)] E ± (0) = 2E f ±8e (14) 



where E±(Q) are the (empirically unknown) zero-CMM energies of the 2e- and 2h-CPs, 
respectively. The quantity Ef is available as a possibly convenient energy scale. It is 
not to be confused with the Fermi energy Ef = \mv F = h 2 k F /2m = k^Tp where Tp 
is the Fermi temperature. If n = N/L d is the total number-density of charge-carrier 
electrons, then n = k F /2n in 2D and = k F /3n 2 in 3D. Thus, the Fermi energy Ep 
equals 7th 2 n/m in 2D and (h 2 /2m)(3K 2 n) 2 ^ in 3D, while Ef equals 7th 2 rif/m in 2D 
and (h 2 /2m)(37t 2 tif) 2 / 3 in 3D, i.e., is the same as Ep with n replaced by n } which will 
serve as a convenient density scale. The quantities Ef and Ep coincide only when perfect 
2e/2h-CP symmetry holds, i.e., when n = rif. 

The interaction Hamiltonian (flOl) can be further simplified by dropping all K 7^ 
terms. This is also done in BCS theory but in the full BCS Hamiltonian (fTT|) . However, in 
the GBEC theory these terms are retained in © so that Hq describes a non-Fermi-liquid 
normal state. Following Bogoliubov |[39|l . the boson operators b^ : bo and Cq , cq remain- 
ing in (flOl) are then replaced by the "c-numbers" ^Nq(T) and ^JMq(T), respectively, 
where Nq and Mo are the number of K = 2e- and 2h-CPs. The full GBEC Hamiltonian 
H thus becomes bilinear in the fermion operators a^ s , a^ s and is diagonalizable via a 
Bogoliubov transformation exactly without assuming weak coupling as in BCS theory. 
One constructs the grand potential Q. for the full GBEC Hamiltonian as 

a(T 1 L d , y i 1 N M = -k B T\n[Tr exp {-j6(# - nN op )}] (15) 

where "7V" stands for "trace" and /3 = XjkpfT- Minimizing with respect to A^o (the 
number of zero-CMM 2e-CPs) and Mq (the same for 2h-CPs), while simultaneously 
fixing the total number N of electrons via the electron chemical potential /1, determines 
an equilibrium state of the system with volume L d and temperature T. One thus requires 
that 

dQ. n dQ. dQ. 

- — =0, - — =0 and —= -N. (16) 

dN dM dpi V ' 

Here evidently includes both paired and unpaired CP electrons. Some algebra then 
leads to the three coupled integral Eqs. (7)-(9) of Ref. fl 1 IN . The relation between the 
fermion spectrum E(e) and fermion energy gap A(e) turns out to be of the BCS form 



E(e) = \ (e-yi) 2 +A 2 (e) butwhere A(e) = v /w^/+(e) + v ^n "/_(e). (17) 



This last expression for the gap A(e) implies a simple ^-dependence rooted in the 2e-CP 
n (T) = N (T) /L d and 2h-CP m (T) = M (T) /L d number densities of BE-condensed 
bosons, i.e., 

A(T) = y/n^T)U{e) + y/ m {T)f-{e). (18) 

This generalizes the relation between BCS and BEC order parameters first reported 
in Ref. |8|]. Self-consistent (at worst, numerical) solution of the aforementioned three 
coupled equations then yields the three thermodynamic variables of the GBEC theory 



no(T,n,[l), mo(r,«,/i) and fl(T,n). 



(19) 



H = H + H inl = H e/h + H e/h p a i rs + { ■ ■ ■ } 
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FIGURE 1. Flowchart outlining conditions under which the GBEC [or "complete boson-fermion 
model" (CBFM)] theory reduces to all five statistical theories of superconductors (ovals). Symbols are 
explained throughout text. 



Most significantly, the three GBEC theory equations contain the key equations of five 
different statistical theories as special cases; for a detailed review see Ref. [14]. Figure 1 
illustrates this in a flowchart. Perfect 2e/2h CP symmetry signifies equal number of 2e- 
and 2h-CPs, i.e., n B {T) = m B (T) as well as n (T) = m (T). With d) Eqs. ([28]) and 
(|29|) below imply that Ef coincides with /i, and the GBEC theory then reduces to: 

i) the gap and number equations of the BCS-Bose crossover picture for the BCS model 
interaction — if the BCS parameters V and Debye energy hcoo are identified with the BF 
interaction Hamiltonian Hi nt parameters f 2 /28e and 8e, respectively. The crossover 
picture for unknowns A(T) and jl{T) is now supplemented by the key expression 
relating BCS with BEC precisely, namely 



A(T)=fyfidT)=fy/mo(T). (20) 

The crossover picture is associated with many authors beginning in 1967 with Friedel 
and coauthors [40] and then given mayor impetus by Eagles 114 1[| who in turn introduced 



the BEC mechanism into the picture; for reviews see Refs. [42] [43]. However, it is 
widely unrecognized to be a very modest improvement, at least for the Cooper/BCS 
model interaction, over BCS theory per se since, e.g., an unphysically large X of about 



8 is required to bring n(T c )/Ep in 2D down from 1.00 to 0.998; indeed, ^-values differ 
very slightly rf23ll in 2D between the crossover and BCS theories all the way up to 
A ~ 50 when the Fermi surface originally pinned at /i disappears by becoming negative 
so that the model interaction breaks down. In fact, room-temperature superconductivity 
is predicted ll23ll by BCS theory in 2D for a material with 7> = 10 3 K, but only for A 
values somewhat larger than 10; these are still too unphysical as they exceed the Migdal 
Il46ll threshold of A > A for ionic-lattice instability. If one imposes that fi(T c ) = Ep 
exactly, as follows from the number equation for weak BF coupling /, the crossover 
picture is well-known to reduce to: 

ii) ordinary BCS theory which is characterized by a single equation, the gap equation 
for any T . Thus, the BCS condensate is precisely a BE condensate whenever both 
u b{T) = mfi(r) and no(T) = mo(T) and the BF coupling / is small. The condensation 
energies of the GBEC and BCS theories coincide exactly for all coupling. 

On the other hand, for no 2h-CPs present the GBEC theory reduces 111 ill also to: 

iii) the BEC BF model in 3D of Friedberg and T.D. Lee |9l llOJ characterized by 
the relation A(T) = fy/no(T) first reported in Ref. Jsl]; but lacking 2h-CPs this model 
cannot be fully related to BCS theory. When f — it reduces to: 



iv) the ideal BF model (IBFM) of Refs. |44j, |45J] that predicts nonzero 2e-CP BEC 



r c 's even in 2D. The "gapless" IBFM cannot describe the superconducting phase. But 
considered as a model for the normal state it should provide feasible T c 's as singularities 



within a BE scenario that are approached from above T c , and this is indeed [45] the case. 
Finally, in the limit of no unpaired electrons this model in 3D reduces to: 

v) the familiar r c -formula of ordinary BEC in 3D but where the boson number-density 
is temperature dependent. 

The vastly more general GBEC theory has been applied in both 2D and 3D and 
gives sizeable enhancements in r c 's over BCS theory for moderate departures from 
perfect 2e/2h-pair symmetry. This is attained for the same Cooper/BCS interaction 
model (coupling strength A and cutoff hcoo) parameter values often used in conventional 
SCs. The three coupled equations of the GBEC theory that determine the J-dimensional 
BE-condensate number-densities no(T) and mo(T) of 2e- and 2h-CPs, respectively, as 
well as the electron chemical potential n(T), were first solved numerically Iil2ll around 
the BCS point of the T c vs. n phase diagram. At n/nf = 1 one has perfect 2e/2h-CP 
symmetry; the plain n/ can be seen to be the number density rtf{T) of unpaired but 
BCS -correlated electrons when A = and T = 0, whenever /i = Ef. In general 

CO 

n f (T) = J deN(e)[\ -^tanh^e)] (21) 



o 



with E(e) defined in (fT7T) . This expression is precisely the BCS expression for the 
electron number density 

n = £v 2 k (T) (22) 



k,.v 



where Vk(T) is the temperature-dependent Bogoliubov function. Alongside two gap- 
like equations involving no(T) and mo(T), the third, or "complete" number, equation 



explicitly reads 



n f (T) +2n (T) +2n B+ (T) -2m (T) -2m B+ (T) = n (23) 

with m B+ (T), e.g., being precisely the number of "pre-formed" K > 2h-CPs, and 
n B +(T) the same for 2e-CPs. Besides the normal phase consisting of the ideal BF 
ternary gas described by Hq, three different stable BEC phases emerge: two pure phases 
consisting of a 2e-CP BEC and a 2h-CP BEC, as well as a mixed phase consisting 
of both types of BECs in varying proportions. For a half-and-half mixed phase all the 
boson number-density terms in (1231 cancel out and the BCS number equation rif(T) = n 
is recovered. 

We shall focus on the linear dispersion that occurs in leading order in K for "ordinary" 
CPs in a Fermi sea as well as for "generalized" CPs in a BCS-correlated state. For the 
latter, the boson excitation energy £ to be used has a leading term in the many-body 
Bethe-Salpeter (BS) CP dispersion relation given by e ~ (X/4)HvfK in 3D Q30Q , as 
part of an expansion similar to © in 2D. As before, A = VN(Ep) where N(Ep) is the 
electron DOS (for one spin) at the Fermi surface. Note that e is no longer the quadratic 
h 2 K 2 /2(2m) often assumed [H - 11211 . 114311 and associated with a composite boson 



of mass 2m moving not in the Fermi sea but in vacuum. 

CRITICAL TEMPERATURES IN 3D: GBEC WITH / = 

Fully equivalent to a "complete boson-fermion model" (CBFM), the GBEC theory with 
/ = becomes an "ideal boson-fermion model" (ICBFM), see Fig. 1. The ICBFM is 
completely described in d dimensions by the Hamiltonian Hq defined by ©. One can 
construct its associated grand potential as 

Clo(T,L d ,n,No,Mo) = -k B T\n[Tr exp{-P(H - ylN op )}] (24) 

with Nq and Mq as before the number of zero-CMM 2e- and 2h-CPs, respectively. One 
gets 

Q (T, h\\i, Nq, Mq) = [E+(0) - 2/i]iV + [2/i -£_(0)]M 

oo 

- 2k B TB d L d J k^dklnll+expi-piEk-ii)}] 
o 

oo 

+ k B TB d L d J K d ~ l dK\n{\ -exp{-/3[£+(^) - 2/i]}) 

o+ 

oo 

+ k B TB d L d J K d ^dK\n(l -exp{-/3[2jii -£_(£)]}) (25) 



o+ 



where B d = \j2% and 1/27T 2 for d = 1, 2 and 3, respectively. An equilibrium 
thermodynamic state makes Q,q(T,L ,H,Nq,Mo) stationary with respect to M) and to 



Mo and requires that the number density of electrons n=N/L d remain constant. Thus, 
one imposes that 

|^=0, |5°=0 and -^=N. (26) 
dN dM d/i 

This leads to the three relations 

[E+(0)-2/j]=0, [2/i-£_(0)] =0 and n f {T) + 2n B (T) -2m B {T) = n (27) 

with the latter being the "number equation" that ensures charge conservation in the 
ternary mixture. In the last relation 

oo 

n B {T) = n (T)+B d J dKK d - l [exp{j5[E + (0)-2^ + e K ]}-\]- 1 (28) 

o+ 

m B (T) ee m (r)+B d y'd^- 1 [exp{jS[2/i- J B_(0) + e z ]-l]- 1 



(29) 



o+ 

while 



n f (T) = 2B d J J^- 1 [exp{/3(e,-Ai)} + l]- 1 . (30) 
o 

This last expression can be interpreted as the number density of unpaired electrons at 
any T and £ k = h 2 k 2 /2m is the electron energy. 

Consider the three equations (1271) assuming only that E (0) < E + (0) . In 3D, the third 
equation of (|27T) explicitly becomes 

oo 

n ee 4/37T 2 = ^ 2 y^[exp{/3(e fc - i u)} + l] _1 
o 

CO 

+2n (T)+7t- 2 J K 2 dK(exp{p[E+(0)-2n + £ K ]}- l)" 1 
o 

CO 

-2m (r)-7zr 2 J K 2 dK(Qxp{P[2[i - E_(0) + e K }} - 1)' 1 . (31) 
o 

Take first the limiting case 2/i = E+(0) when the first equation of (|27T) is satisfied 
but not the second. From (l26l) this implies that Mq(T) must vanish for all T, but that 
Nq(T) 7^ at least for some 7 1 so that a pure 2e-CP BEC phase may occur below a critical 
temperature T c (possibly zero) determined by no(T Cl n) = 0. After substituting (fl4l) in 
(ED with ju = £ + (0) /2, thus eliminating £ ± (0) in favor of E f = (h 2 /2m) (3n 2 n f ) 2 / 3 and 
8e ee we obtain the dimensionless "working number equation" for the pure 2e-CP 



BEC phase critical temperature 



1/3 = Jk 2 dk[exp{(k 2 -n- 2 / 3 -v/2)/T c } + l}~ 1 + J K 2 dK[exp(XK/2T c ) - l]' 1 
o o 

oo 

-J K 2 dK[Qxp{{XK + Av)/2f c }-l} 1 (32) 



o 

where k = kjkp^K^ = K/kp, h = n/tif, T c = T c /Tp, V = hCQo/Ep, and we took ~ 
XhvpK/A [Ref. ri30h - Eq. (12)] for d = 3. The integrals are exact, the first and last being 
expressible as polylogarithm functions Li a (z) or PolyLog[a, z] rf47ll where 

-aL, CT (-«z) = ^-r^ = —aL— (33) 

with z an effective fugacity. For a = —I (|33l) reduces to the Bose integral gcr(-z) which 
for z = 1 and c > 1 becomes the Riemann Zeta function £(cr) of order a; for a = 1 
(1331 becomes the Fermi integral /<j(z). Both integrals are as defined in Appendices D 
and E of Ref. [48] . Since Li a {\) = ^(ct), the second integral in (|32l) gives C(3), and the 
working number equation simplifies to 

1/3 = - ^ PoZyLog[3/2, - exp{(v/2 + n- 2 / 3 )/T c }] 
\6T 3 

+ -j§-{a3)-PolyLog[3,exp(-2v/T c )]}. (34) 

This can now be solved for t c as a function of n which is plotted as the dashed curve in 
Fig. 2 for A = \ and v = 0.005. 

The pwre 2h-CP BEC phase comes from the limiting case 2/1 = E-(0) as then the 
second, but not the first, relation in (1271) is satisfied, so that, again from (1261 ), Nq(T) = 
for all r but Mq(T) ^ for some T. A working number equation similar to but different 
from (T3~2l) follows and the critical temperature for this phase is now determined by 
mo(T c ,n) = 0. It eventually gives 

rat 3/2 

1/3 = -^^Po/3;Lo^[3/2,-exp{(n- 2 / 3 -v/2)/f c .}] 
16? 3 

+ -j^-{PolyLog[3,exp(-2v/T c )}-C(3)}- (35) 

No bounded solution of this equation for T c when X = 1/2 and V = h(Oo/E F = 0.005 
was found. The BCS value from their formula T c /Tp ~ l.l34(ha>D/Ef )exp( — 1/A) is 
~ 0.0008 (large black dot in Fig. 2); it lies within the range T c /Tp « 10 -3 empirically 
found for conventional, elemental SCs. Empirical data Q49Q for "exotic" SCs, however, 
fall within the range T c /Tp ~ 0.01 —0.05. Thus, moderate departures from perfect 2e/2h- 
CP symmetry enable the ICBFM to access, unlike BCS theory, empirical T c values for 
exotic SCs without abandoning electron-phonon dynamics. 




= J conventional SCs 
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FIGURE 2. Temperature T (in units of 7>) vs. electron density n (in units of n/ as defined in text) phase 
diagram in 3D showing the critical temperature T c phase boundary for the pure 2e-CP BEC phase (dashed 
curve) for A = 1 /2 with hco D /E F = 0.005. The BCS T C /T F ~ 1.134(hco D /E F )e\p(-l/X) gives ~ 0.0008 
for these values of X and hcOo/Ep; it is marked by a large dot. There is no pure 2h-CP BEC phase solution 
for njnf > 1, nor a mixed phase with both types of BECs for any n/rif > 0, as reported in Ref. Ill 211 for 
f 7^0 but for a quadratic CP dispersion. RTSC refers to room-temperature superconductivity in a material 
with T F = 10 3 K. 



Finally, for intermediate values of /!, namely for E-(0) < 2ji < E+(0), neither the 
first nor second equations of (|26|) are satisfied so that hq(T) = = mo(T) for all T; this 
implies no condensed phases whatsoever. Thus, the ICBFM (characterized by zero BF 
coupling /) contains no mixed phase in contrast to the CBFM where / ^ (Ref. [fl2ll 
for a quadratic dispersion). This case will be treated elsewhere with the correct linear 
dispersion. 

CONCLUSIONS 

Cooper pairs (CPs) are meaningless if referred to the ideal Fermi gas "sea" when 
hole pairs are included along with electron pairs, but survive as positive-energy, finite- 



lifetime, plasmon-like objects with a linear (instead of quadratic) rise in total, or center- 
of-mass, momentum K when referred to a BCS-correlated sea instead. 

The new generalized BEC (GBEC) theory includes as limiting cases the following 
theories: i) BCS and ii) BCS-Bose "crossover," when the BE condensate consists of 
equal numbers of electron- and hole-pairs. It also contains: iii) the Friedberg-T.D. 
Lee BEC model, iv) the "ideal boson-fermion model" (IBFM), and v) ordinary BEC 
theory when there are no unpaired fermions. The BCS condensate is precisely a BE 
condensate of equal numbers of 2e/2h-pairs and weak coupling. Without abandoning 
electron-phonon dynamics the GBEC theory leads to 2-to-3 order-of-magnitude higher 
T c 's — including room-temperature superconductivity. All this rests on four essential 
ingredients: 1) 2h-CPs cannot and must not be neglected in a fully self-consistent 
treatment in any many-fermion system, otherwise a spurious value of T c may result that 
corresponds not to a stable but rather to a metastable state; 2) CPs are bosons, even 
though BCS pairs not; 3) CPs are linearly-dispersive for small K; 4) to achieve higher 
T c 's one must depart from the perfect 2e-/2h-CP symmetry of the BCS condensate which 
in fact is a BE condensate. Neglected in the GBEC theory thus far, however, are: a) 
K > terms in the boson-fermion vertex interactions; b) boson-boson interactions (as 
also in BCS theory); c) a T > Bethe-Salpeter CP treatment; d) different hole and 
electron effective masses; and e) ionic-lattice crystallinity effects which might initially 
be introduced via Van Hove singularities in the electronic DOS or via "bipolarons" 
instead of CPs. 

Finally, at least two mysteries have surfaced here. In spite of neglecting boson- 
boson interactions between severely-overlapping CPs, why has the BCS theory been 
so successful in describing at least conventional SCs? Why are simple models (such as 
the BCS or the GBEC theories) quite able to be of any relevance whatsoever in such 
complex strongly-interacting many-electron systems like SCs? 
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